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CONTINUUM  SHAPE  SENSITIVITY  FOR 
NONLINEAR  AEROELASTIC  GUST  RESPONSE 

Robert  A.  Canfield 


Abstract 

The  objective  of  this  research  was  to  derive  and  apply  new  computational  algorithms  for  the 
aeroelastic  design  sensitivity  of  high  fidelity,  coupled  fluid-structure  interaction  (FSI)  problems. 
A  continuum  sensitivity  analysis  (CSA)  method  was  developed  for  efficient  and  accurate 
computation  of  shape  design  derivatives  for  nonlinear  aeroelastic  structures.  The  local  and  total 
forms  for  CSA  were  compared  for  the  first  time,  based  on  their  accuracy,  efficiency  and 
implementation.  A  proof  was  published  for  equivalence  of  the  total-form  CSA  method  and  the 
discrete  analytic  method — for  the  first  time  for  general  nonlinear  transient  second-order  partial 
differential  equations  that  govern  many  field  problems.  Accuracy  of  the  local-form  CSA  was 
increased  in  two  ways.  One  was  to  use  higher-order  p-elements.  The  other  way  was  to  use  spatial 
gradient  reconstruction  (SGR)  of  the  analysis  output  to  compute  the  sensitivity  boundary 
conditions.  SGR  for  the  structural  domain  was  improved  by  reconstructing  gradients  of  the 
Neumann  sensitivity  boundary  conditions  directly  from  force  and  moment  output  of  the 
structural  analysis  instead  of  reconstructing  higher-order  derivatives.  Guidelines  were  developed 
for  defining  SGR  patches  with  controllable  accuracy  and  rate  of  convergence  for  CSA. 

Summary 

AFOSR  awarded  the  grant  in  June  2014.  This  grant  builds  upon  research  begun  in  AFOSR  grant 
begun  in  FY 1 1 .  The  three-year  follow-on  research  effort  has  further  developed  the  computational 
algorithms  for  efficient  analysis  and  design  optimization  of  aeroelastic  models  in  several  areas. 
Domain  ad  boundary  velocity  approaches,  also  known  as  total  and  local  derivative  formulations, 
were  compared  thoroughly.  The  remedy  for  handling  interface  discontinuities  in  built-up 
structures  for  the  boundary  velocity  CSA  approach  was  extended  from  beams  to  unstiffened  and 
stiffened  plates.  Spatial  gradient  reconstruction  (SGR)  was  demonstrated  as  a  means  to  improve 
CSA  accuracy  in  a  controllable  and  nonintrusive  manner.  Three  journal  articles  [1-3]  were 
published  and  four  others  [4-7]  submitted  for  publication  documenting  the  funded  research. 
Seven  conference  papers  were  presented  during  the  performance  period  [8-14].  Two  dissertations 
were  published:  http://hdl.handle.net/ 1091 9/23282  and  http://hdl.handle.net/10919/48589. 

Progress  was  made  in  the  first  year  in  developing  continuous  sensitivity  equations  (CSE)  for 
fluid-structure  interaction  (FSI).  The  CSE  and  their  associated  sensitivity  boundary  conditions 
were  studied  for  a  built-up  joined  beam  structure  under  transient  aerodynamic  loads.  Eor  built-up 
structures,  the  state  variable  field  and  the  local  sensitivity  variable  field  have  different  continuity 
requirements — the  former  is  C*’  continuous,  while  the  latter  is  piecewise  C*’  continuous.  Eor  the 
discrete  system,  this  might  imply  a  choice  of  different  degrees  of  freedom,  which  in  turn  would 
alter  the  assembly  procedure  for  the  system  matrices.  A  remedy  was  provided  so  that  for 
nonlinear  problems,  when  the  Newton-Raphson  method  is  used,  the  tangent  stiffness  matrix  may 
still  be  used  as  the  sensitivity  coefficient  matrix  for  solving  the  linear  sensitivity  equations  with 
the  Galerkin  finite  element  method.  The  domain  velocity  (total  derivative)  formulation  of  CSE, 
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most  often  used  for  structural  problems  in  part  to  avoid  local  sensitivity  discontinuity  due  to 
strain  discontinuities,  and  the  boundary  velocity  (local  derivative)  formulation,  most  often  used 
for  fluid  problems,  were  compared.  The  discretized  domain  velocity  formulation  of  the  CSE 
method  was  shown  to  be  equivalent  to  the  discrete  analytical  sensitivity  formulation  for  potential 
flow  problems  for  the  first  time.  As  a  result,  the  boundary  velocity  formulation  for  CSE  is  being 
pursued  for  ESI  problems.  Benchmark  examples  with  analytical  and  finite  difference  sensitivity 
results  were  solved  with  CSE  for  the  purpose  of  verification.  The  boundary  velocity  CSE 
formulation  was  implemented  and  solved  in  a  “black  box”  solver  without  access  to  the  source 
code  in  order  to  numerically  solve  the  sensitivity  systems  of  a  nonlinear  ESI  problem  involving 
determinate  and  indeterminate  (built-up)  structures.  Sources  of  inaccuracy  in  the  boundary 
velocity  formulation  relative  to  the  domain  velocity  formulation  were  investigated. 

In  the  second  year  progress  was  made  in  developing  a  continuum  sensitivity  method  for  efficient 
and  accurate  computation  of  shape  design  derivatives  for  nonlinear  aeroelastic  structures  subject 
to  transient  aerodynamic  loads.  The  CSEs  are  often  posed  in  terms  of  local  derivatives  (local 
form)  for  fluid  applications  and  in  terms  of  total  derivatives  (total  form)  for  structural 
applications.  The  local  form  and  the  total  forms  for  static  and  transient  aeroelastic  (fluid- 
structure  interaction)  problems  were  compared  for  the  first  time,  based  on  their  accuracy, 
efficiency  and  implementation.  The  total  form  CSE  exhibited  higher  accuracy  than  local  form, 
when  the  same  numerical  method  was  used  for  solving  fluid  and  structures  benchmark  problems. 
The  local  form  requires  design  velocity  only  at  the  boundaries,  and  the  continuity  requirements 
for  the  design  velocity  is  lower  than  the  total  form  CSE  method.  On  the  other  hand,  local  form 
requires  higher  order  spatial  derivatives  to  specify  its  sensitivity  boundary  conditions  than  does 
total  form  CSE.  Accuracy  of  the  local  form  CSE  was  increased  in  two  ways.  One  was  to  use 
higher-order  /»-elements.  The  other  way  was  to  use  spatial  gradient  reconstruction  of  the  analysis 
output  to  compute  the  sensitivity  boundary  conditions.  The  latter  was  improved  for  the  structural 
domain  by  reconstructing  gradients  of  the  Neumann  sensitivity  boundary  conditions  directly 
from  force  and  moment  output  of  the  structural  analysis  instead  of  reconstructing  higher-order 
derivatives  of  displacement.  The  total  form  continuum  sensitivity  method  was  applied  for  the 
first  time  to  nonlinear  transient  gust  response.  The  equivalence  of  the  discrete  analytic  method 
and  the  total  form  CSE  method  was  proved  for  the  first  time  for  general  nonlinear  transient 
second-order  partial  differential  equations  that  govern  many  field  problems.  The  equivalence  of 
these  two  methods  was  demonstrated  for  general  transient  geometric  nonlinear  elasticity  and 
potential  flow,  as  specific  examples. 

In  the  third  year,  continuum  sensitivity  analysis  (CSA)  was  matured  for  efficient  and  accurate 
computation  of  shape  design  derivatives  for  nonlinear  aeroelastic  structures  subject  to  transient 
aerodynamic  loads.  An  advantage  of  the  local  form  is  that  it  requires  design  velocity  only  at  the 
boundaries,  and  the  continuity  requirements  for  the  design  velocity  is  lower  than  the  total-form 
CSA  method.  On  the  other  hand,  local  form  requires  higher  order  spatial  derivatives  to  specify 
its  sensitivity  boundary  conditions  than  does  total-form  CSA.  Accuracy  improvement  of  the 
local-form  CSA  was  increased  by  way  spatial  gradient  reconstruction  (SGR)  of  the  analysis 
output  to  compute  the  sensitivity  boundary  conditions  was  published  for  structures  for  the  first 
time.  An  advance  compared  to  the  previous  application  for  fluids,  SGR  for  the  structural  domain 
was  improved  by  reconstructing  gradients  of  the  Neumann  sensitivity  boundary  conditions 
directly  from  force  and  moment  output  of  the  structural  analysis  instead  of  reconstructing  higher- 
order  derivatives  of  displacement.  Studies  this  year  indicated  that  CSA  with  SGR  achieves  the 
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same  rate  of  convergence  as  the  original  analysis,  provided  sufficient  layers  are  employed  in  the 
patches  for  Taylor  series  reconstruction.  An  example  demonstrated  that  local-form  CSA  with 
SGR  can  approach  the  accuracy  achieved  using  analytical  spatial  derivatives  and  can  exceed  that 
accuracy  of  discrete  sensitivity  or  equivalently  total-form  CSA.  Guidelines  were  developed  for 
defining  SGR  patches  with  controllable  accuracy  and  rate  of  convergence. 

Accomplishments 

Several  significant  contributions  to  the  state  of  knowledge  were  made  for  the  application  of 
CSA.  Although  they  were  posed  in  support  of  the  sensitivity  analysis  of  FSI  problems,  some 
findings  are  more  general  and  more  fundamental  than  the  application  to  FSI. 

The  total  form  continuum  sensitivity  method  was  applied  for  the  first  time  to  nonlinear  transient 
aeroelastic  problems.  For  nonlinear  problems,  it  was  shown  that  the  tangent  stiffness  matrix 
yields  the  desired  sensitivity  coefficient  matrix  for  solving  the  sensitivity  equations.  The  local 
form  CSE  was  applied  for  the  first  time  to  built-up  structures.  The  difficulty  in  implementation 
was  pointed  out  and  resolved  [1].  The  first  comparison  of  the  local  form  and  the  total  form  CSE 
method  for  static  and  transient  aeroelastic  problems  has  been  made  on  their  accuracy,  efficiency 
and  implementation  [6].  The  equivalence  of  the  discrete  analytic  method  and  the  total  form  CSE 
method  is  proved  for  the  first  time  for  general  nonlinear  transient  second-order  partial 
differential  equations  that  govern  many  field  problems  [2]. 

The  total  form  CSE  method,  once  discretized,  was  proven  to  be  equivalent  to  the  discrete 
analytic  method  for  nonlinear  second-order  field  problems  under  certain  conditions.  The 
conditions  are  that  the  two  methods  use  the  same:  (1)  numerical  approximation  of  the  field 
variables  and  sensitivity  variables,  (2)  integration  of  the  finite  element  matrices,  (3)  design 
velocity  fields  that  are  linear  with  respect  to  the  design  variable,  and  (4)  shape  functions  for 
domain  transformation  and  for  design  velocity  field  calculation.  Although,  under  these 
conditions,  the  discrete  analytic  method  is  identical  to  the  total  form  CSE  method,  the  total  form 
CSE  can  be  implemented  without  satisfying  all  these  conditions.  The  advantage  of  the  CSE 
method  is  that  the  discrete  analytic  method  requires  the  knowledge  of  the  implementation  details 
of  the  original  analysis  codes,  which  may  be  inaccessible  or  very  complex,  while  the  continuum 
method  requires  minimum  information  of  the  source  codes.  The  continuum  method  only  requires 
knowledge  of  the  governing  equations,  the  numerical  solutions  and  the  spatial  gradients  of  the 
solutions  for  the  governing  equations.  Thus,  the  CSE  method  can  be  implemented  by  treating  the 
existing  analysis  solver  as  a  “black  box”  [3]. 

A  primary  objective  of  this  research  was  to  develop  design  sensitivity  analysis  (DSA)  that  can 
accurately  and  efficiently  calculate  design  derivatives  of  high  fidelity  models,  such  as  nonlinear 
transient  aeroelastic  gust  analysis.  Eurthermore,  the  DSA  method  should  be  nonintrusive, 
providing  the  ability  to  calculate  design  derivatives  of  analyses  conducted  using  general-purpose 
codes.  This  specific  DSA  formulation  developed  here  is  called  local  CSA  with  spatial  gradient 
reconstruction  (SGR).  Moreover,  an  added  benefit  of  using  local  CSA  is  that  mesh  sensitivities 
are  not  required.  Specific  contributions  follow. 

•  This  research  marks  the  first  application  of  local  CSA  for  linear  [3]  and  nonlinear  analysis  [4] 
of  structural  models  using  plate  elements.  It  is  the  first  demonstration  of  implementing  local 
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CSA  with  black  box  analysis  tools  without  knowledge  of  finite  element  shape  funetions. 
Furthermore,  it  is  the  first  application  of  local  CSA  to  a  built-up  structure  that  uses  mixed 
element  types  (plate  and  beam  elements).  Interface  boundary  conditions  were  derived  for 
joined  mixed  elements. 

•  Local  CSA  with  SGR  can  be  implemented  with  general  purpose  codes  in  a  nonintrusive 
fashion,  treating  the  analysis  tool  as  a  “black  box.”  The  only  requirements  are  for  the  user  to 
differentiate  the  boundary  conditions  to  define  sensitivity  loads  and  for  the  analysis  tool  to 
provide  as  output  the  system  response  of  primary  and  seeondary  variables,  as  well  as  the 
system  matriees.  The  method  is  element  agnostie  and  can  be  readily  implemented  for  a  variety 
of  analysis  models  with  varying  degrees  of  fidelity.  Numerous  examples,  conducted  with 
Nastran  analyses,  demonstrated  these  attributes.  The  results  also  demonstrate  that  the  method 
can  be  as  aeeurate  as  other  analytic  DSA  methods  and  more  aeeurate  than  most  [5]. 

•  SGR  has  been  utilized  in  the  context  of  flow  sensitivities,  but  this  is  the  first  use  of  SGR  in  the 
context  of  structural  sensitivities.  Previously,  SGR  was  utilized  to  approximate  high-order 
derivatives  of  primary  variables  to  formulate  local  CSE  boundary  conditions.  The  local  CSA 
method  with  SGR  presented  here  requires  only  the  approximation  of  first-order  derivatives,  but 
of  both  primary  and  seeondary  variables.  This  is  a  key  ingredient  in  making  the  method 
element-agnostic  and  amenable  to  general  purpose  codes.  This  research  marks  the  first  time 
that  SGR  has  been  employed  in  this  way  and  for  this  purpose  [3]. 

•  Guidelines  were  developed  to  control  convergence  of  aeouraey  for  loeal  CSA  with  spatial 
gradient  reeonstruetion.  In  general,  the  eonvergenee  rate  is  largely  dependent  on  the  aeeuraey 
of  the  spatial  derivative  approximations  that  appear  in  the  local  CSE  boundary  conditions.  The 
accuracy  of  these  approximations  is  driven  by  the  parameters  used  for  the  SGR  [5]. 

•  A  hybrid  adjoint  was  derived  to  effieiently  compute  statie  shape  derivatives  for  many  design 
variables,  using  a  diserete  adjoint  with  CSA  for  whieh  aeouraey  is  oontrolled  by  SGR  [7,  14]. 

Acknowledgment/Disclaimer 

This  work  was  sponsored  by  the  Air  Eorce  Office  of  Scientific  Research,  USAE,  under 
grant/oontract  number  FA9550-1 1-1-01 17.  The  views  and  conclusions  contained  herein  are 
those  of  the  authors  and  should  not  be  interpreted  as  neoessarily  representing  the  offioial  polioies 
or  endorsements,  either  expressed  or  implied,  of  the  Air  Force  Office  of  Scientific  Research  or 
the  U.S.  Government. 

References 

[1]  Liu,  Shaobin;  Canfield,  Robert.  A.  (2013),  “Boundary  Veloeity  Method  for  Continuum  Shape 
Sensitivity  of  Nonlinear  Fluid-Structure  Interaction  Problems,”  Journal  of  Fluids  and  Structures, 
http://dx.doi.Org/10.1016/j.jfluidstructs.2013.05.003,  June  6,  2013,  Vol.  40,  pp.  284-301 

[2]  Liu,  Shaobin;  and  Canfield,  Robert.  A.  (2013),  “Equivalenee  of  Continuum  and  Discrete 
Analytic  Sensitivity  Methods  for  Nonlinear  Differential  Equations,”  Struetural  and 
Multidisciplinary  Optimization,  doi:  10.1007/s00158-013-0951-4  (online  4  July  2013), 

Deeember  2013,  Volume  48,  Issue  6,  pp.  II73-II88 


AEOSR  GRANT  #FA9550-I  I-I-OI 17 


5 


[3]  Cross,  David  M.  and  Canfield,  Robert  A.  (2014),  “Loeal  eontinuum  shape  sensitivity  with 
spatial  gradient  reeonstruetion,”  Structural  and  Multidisciplinary  Optimization,  June,  2014, 
DOI  10.1007/S00158-014-1092-0 


[4]  Cross,  David  M.;  and  Canfield,  Robert  A.,  “Loeal  Continuum  Shape  Sensitivity  with  Spatial 
Gradient  Reeonstruetion  for  Nonlinear  Analysis,”  aeeepted  for  publieation  in  Structural  and 
Multidisciplinary  Optimization 

[5]  Cross,  David  M.;  and  Canfield,  Robert  A.,  “Convergenee  of  Aeeuracy  for  Loeal  Continuum 
Sensitivity  Method  Using  Spatial  Gradient  Reeonstruetion,”  submitted  to  AIAA  Journal 

[6]  Liu,  Shaobin;  and  Canfield,  Robert  A.,  “Two  Forms  of  Continuum  Shape  Sensitivity  Method 
for  Fluid-Strueture  Interaetion  Problems,”  submitted  to  AIAA  Journal 

[7]  Kulkami,  Mandar  D.;  Cross,  David  M.;  and  Canfield,  Robert  A.,  “Diserete  Adjoint 
Formulation  for  Continuum  Sensitivity  Analysis,”  submitted  to  AIAA  Journal 

[8]  Liu,  Shaobin;  and  Canfield,  Robert  A.,  “Continuum  Shape  Sensitivity  Method  for  Fluid  Flow 
Around  an  Airfoil,”  53rd  AIAA/ASME/ASCE/AHS/ASC  Struetures,  Struetural  Dynamics  and 
Materials  Conference,  AIAA-2012-1426,  Honolulu,  Hawaii,  Apr.  23-26,  2012 

[9]  Cross,  David;  and  Canfield,  Robert  A.,  “Solving  Continuum  Shape  Sensitivity  with  Existing 
Tools  for  Nonlinear  Aeroelastic  Gust  Analysis,”  53rd  AIAA/ASME/ASCE/AHS/ASC 
Structures,  Structural  Dynamics  and  Materials  Conference,  AIAA-2012-1923,  Honolulu, 

Hawaii,  Apr.  23-26,  2012 

[10]  Eiu,  Shaobin;  and  Canfield,  Robert  A.,  “Continuum  Sensitivity  Method  for  Aeroelastic 
Shape  Design  Problems,”  AIAA-20 12-5480,  14th  AIAA/ISSMO  Multidisciplinary  Analysis  and 
Optimization  Conference,  September  17-19,  2012 

[11]  Cross,  David;  and  Robert  Canfield,  “Continuum  Shape  Sensitivity  with  Spatial  Gradient 
Reconstruction  of  Nonlinear  Aeroelastic  Gust  Response,”  AIAA-2012-5597,  14th 
AIAA/ISSMO  Multidisciplinary  Analysis  and  Optimization  Conference.  September  17-19,  2012 

[12]  Cross,  David  M.;  and  Canfield,  Robert  A.,  “Continuum  Shape  Sensitivity  with  Spatial 
Gradient  Reconstruction  of  Built-up  Structures,”  DOI:  10.25 14/6. 2013-1933,  54th 
AIAA/ASME/ASCE/AHS/ASC  Structures,  Structural  Dynamics,  and  Materials  Conference, 
April  8-11,2013 

[13]  Kulkarni,  Mandar  D.;  Canfield,  Robert  A.;  and  Path,  Mayuresh,  “Integration  of  Geometric 
Sensitivity  and  Spatial  Gradient  Reconstruction  for  Aeroelastic  Shape  Optimization,”  AIAA- 
2014-0470,  AIAA  SciTech  lOth  Multidisciplinary  Design  Optimization  Conference,  National 
Harbor,  Maryland,  January  13-17,  2014 

[14]  Kulkarni,  Mandar  D.;  Canfield,  Robert  A.;  and  Path,  Mayuresh,  “Nonintrusive  Continuum 
Sensitivity  Analysis  for  Aerodynamic  Shape  Optimization,”  AIAA-20 14-2043,  15th 
AIAA/ISSMO  Multidisciplinary  Analysis  and  Optimization  Conference,  Atlanta,  Georgia, 

June  16,  2014 


AEOSR  GRANT  #PA9550-1 1-1-0117 


6 


